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Abstract 

Qh After the successful determination of the reactor neutrino mixing angle 

r-| ^13 = 0.16 7^ 0, a new feature suggested by the current neutrino oscillation 

data is a sizeable deviation of the atmospheric neutrino mixing angle ^23 
from 7r/4. Using the fact that the neutrino mixing matrix U = UeUu, where 
>■ Ue and Ui, result from the diagonalisation of the charged lepton and neu- 

^ trino mass matrices, and assuming that Ui, has a i) bimaximal (BM), ii) 

tri-bimaximal (TBM) form, or else iii) corresponds to the conservation of 
O the lepton charge L' = Lg — — L^- (LC), we investigate quantitatively 

04 what are the minimal forms of Ue, in terms of angles and phases it contains, 

that can provide the requisite corrections to Ui, so that ^13, 623 and the so- 
lar neutrino mixing angle 612 have values compatible with the current data. 
Two possible orderings of the 12 and the 23 rotations in Ue, "standard" 
and "inverse", are considered. The results we obtain depend strongly on 
the type of ordering. In the case of "standard" ordering, in particular, the 
Dirac CP violation phase 5, present in U, is predicted to have a value in a 
narrow interval around i) (5 = vr in the BM (or LC) case, ii) 6 = 37r/2 or ir/2 
in the TBM case, the CP conserving values (5 = 0, vr, 2n being excluded in 
the TBM case at more than 4a. 
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1 Introduction 



Understanding the origin of the patterns of neutrino masses and mixing, emerging 
from the neutrino oscillation, /3— decay, etc. data is one of the most challenging 
problems in neutrino physics. It is part of the more general fundamental problem 
in particle physics of understanding the origins of flavour, i.e., of the patterns of 
the quark, charged lepton and neutrino masses and of the quark and lepton mixing. 

At present we have compelling evidence for the existence of mixing of three light 
massive neutrinos z/j, z = 1, 2, 3, in the weak charged lepton current (see, e.g., [1]). 
The masses mj of the three light neutrinos Vi do not exceed approximately 1 eV, 
rrii^l eV, i.e., they are much smaller than the masses of the charged leptons and 
quarks. The three light neutrino mixing we will concentrate on in the present ar- 
ticle, is described (to a good approximation) by the Pontecorvo, Maki, Nakagawa, 
Sakata (PMNS) 3x3 unitary mixing matrix, f/pMNS- In the widely used standard 
parametrisation [1], t/pMNS is expressed in terms of the solar, atmospheric and 
reactor neutrino mixing angles 6'i2, and 6'i3, respectively, and one Dirac - 5, 
and two Majorana [2] - 0^21 and 0:31, CP violation (CPV) phases: 



t^PMNS = U = V{di2, O23, 6*13, <5(a21, "31; 



where 



1 \ / ci3 si^e-'' \ ( C12 S12 
= I C23 S23 1 -S12 C12 I , (2) 



-S23 C23 / V -si^e'' ci3 / V 1 

Q = diag(l, e""''/^, e'"^^/2) _ ^3) 

and we have used the standard notation Cij = cos 9ij, Sij = sin 9ij with < 9ij < 
7'"/2, < (5 < 27r and, in the case of interest for our analysis, < aji < 2tt, j = 2, 3 
(see, however, [3]). If CP invariance holds, we have 5 = 0, tt, and [4] 021(31) = 0, vr. 

The neutrino oscillation data, accumulated over many years, allowed to deter- 
mine the parameters which drive the solar and atmospheric neutrino oscillations, 
Am^i, 6*12 and jAmg^^l = |Am32|, 6*23, with a high precision (see, e.g., [5]). 

Furthermore, there were spectacular developments in the last 1.5 years in what 
concerns the angle ^13 (see, e.g., [1]). They culminated in a high precision deter- 
mination of sin^ 26*13 in the Daya Bay experiment with reactor i/g [6] : 



sm 



2 2^13 = 0.089 ±0.010 ±0.005. (4) 



Similarly the RENO, Double Chooz, and T2K experiments reported, respectively, 
4.9(7, 2.9(7 and 3.2a evidences for a non-zero value of 6*13 [7], compatible with the 
Daya Bay result. The high precision measurement on ^13 described above and 
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Parameter 



Best fit la range 



2a range 



3(7 range 



(NH) 
(IH) 



0.155 0.147- 0.163 

0.307 0.291 - 0.325 

0.386 0.365 - 0.410 

0.392 0.370 - 0.431 

3.39 2.42 - 4.27 

3.42 2.61 - 4.62 



0.353 



0.139 - 0.170 
0.275 - 0.342 
0.348 - 0.448 
0.484 © 0.543 



0.641 



0.130 - 0.177 
0.259 - 0.359 
0.331 - 0.637 
0.335 - 0.663 



Table 1: Summary of the resufts of the global fit of the PMNS mixing angles 
taken from [8] and used in our analysis. The results on the atmospheric neutrino 
angle 623 and on the Dirac CPV phase 5 depend on the type of neutrino mass 
hierarchy. The values of sin^ 623 and 6 obtained in both the cases of normal 
hierarchy (NH) and inverted hierarchy (IH) are shown. 



the fact that ^13 turned out to have a relatively large value, have far reaching 
implications for the program of research in neutrino physics (see, e.g., [1]). After 
the successful measurement of ^13, the determination of the absolute neutrino 
mass scale, of the type of the neutrino mass spectrum, of the nature - Dirac or 
Majorana, of massive neutrinos, as well as getting information about the status of 
CP violation in the lepton sector, are the most pressing and challenging problems 
and the highest priority goals of the research in the field of neutrino physics. 

A global analysis of the latest neutrino oscillation data presented at the Neu- 
trino 2012 International Conference [5], was performed in [8]. The results on 
sin^^i2, sin^ 6'23 and sin^ 6'i3 obtained in [8], which play important role in our fur- 
ther discussion, are given in Table 1. An inspection of Table 1 shows that, in 
addition to the nonzero value of 6*13, the new feature which seems to be suggested 
by the current global neutrino oscillation data is a sizeable deviation of the angle 
623 from the value 7r/4. This trend is confirmed by the results of the subsequent 
analysis of the global neutrino oscillation data performed in [9]. 

Although ^13 ^ 0, $23 ^ tt/A and 612 ^ tt/4, the deviations from these values are 
small, in fact we have sin 6*13 ^ 0.16 < 1, 7r/4 - 6*23 = 0.11 and 7r/4 - 6*12 = 0.20, 
where we have used the relevant best fit values in Table 1. The value of ^13 and 
the magnitude of deviations of 6*23 and 6^12 from 7r/4 suggest that the observed 
values of 613, 623 and 612 might originate from certain "symmetry" values which 
undergo relatively small (perturbative) corrections as a result of the corresponding 
symmetry breaking. This idea was and continues to be widely explored in attempts 
to understand the pattern of mixing in the lepton sector (see, e.g., [10, 11, 12, 13, 
14, 15, 16]). Given the fact that the PMNS matrix is a product of two unitary 
matrices, 

U = UlU,, (5) 
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where Ue and f/^, result respectively from the diagonalisation of the charged lepton 
and neutrino mass matrices, it is usually assumed that Uu has a specific form 
dictated by a symmetry which fixes the values of the three mixing angles in Ui, 
that would differ, in general, by perturbative corrections from those measured in 
the PMNS matrix, while Ue (and symmetry breaking effects that we assume to 
be subleading) provide the requisite corrections. A variety symmetry forms of Ut, 
have been explored in the literature on the subject (see, e.g., [17]). In the present 
study we will consider three widely used forms, 
i) Tribimaximal Mixing (TBM) [18]: 



TBM 



V 







\ 



(6) 



ii) Bimaximal Mixing (BM) [19]: 



BM 



v 



1 

V2 

_ 1 
2 

1 
2 
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V2 
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(7) 



iii) the form of U^, resulting from the conservation of the lepton charge L' 
Lf. — L^ — Lj. of the neutrino Majorana mass matrix [20] (LC): 



where C23 



LC 



cos 6^23 and 



_ ^ 

^23 
V2 



1 

V2 





\ 


''23 






V2 


■=■23 




*23 


•-23 




V2 


/ 



sin ^23- 



We will define the assumptions we make on f/g and in full generality in Sec- 
tion 2. Those assumptions allow us to cover, in particular, the case of corrections 
from Ue to the three widely used forms of U^, indicated above. We would like to 
notice here that if ?7e = 1, 1 being the unity 3x3 matrix, we have: 



in all three cases of interest of U,, 



IJ 1713 

ii) ^23 = 7r/4, if U^ coincides with ?7tbm or U-qmi while 6*23 can have an arbitrary 
value if Uy is given by [/lc! 

iii) 9i2 = 7r/4, for U,y = Ubm or [/lc, while 6*12 = sin"^(l/A/3) if U,y = Utbm- 
Thus, the matrix Ue has to generate corrections 
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i) leading to ^13 7^ compatible with the observations in all three cases of 1/^ 
considered; 

ii) leading to the observed deviation of 6*12 from 7r/4 in the cases of Uu = Ubm or 
Ulc. 

iii) leading to the sizeable deviation of 623 from n / 4 ioi U,, = U^bm ot Ubm, if it 
is confirmed by further data that sin^ 623 = 0.40. 

In the present article we investigate quantitatively what are the "minimal" 
forms of the matrix Ug in terms of the number of angles and phases it contains, 
that can provide the requisite corrections to Utbm, Ubm and [/lc so that the angles 
in the resulting PMNS matrix have values which are compatible with those derived 
from the current global neutrino oscillation data, Table 1. Our work is a natural 
continuation of the study some of us have done in [15] and earlier in [11, 12, 13, 14]. 

The paper is organised as follows. In Section 2 we describe the general setup 
and we introduce the two types of "minimal" charged lepton "rotation" matrix 
Ue we will consider: with "standard" and "inverse" ordering. The two differ by 
the order in which the 12 and 23 rotations appear in Ug. In the same Section 
we derive analytic expressions for the mixing angles and the Dirac phase 6 of the 
PMNS matrix in terms of the parameters of the charged lepton matrix Ug both 
for the tri-bimaximal and bimaximal (or LC) forms of the neutrino "rotation" 
matrix U^- In Sections 3 and 4 we perform a numerical analysis and derive, in 
particular, the intervals of allowed values at a given CL. of the neutrino mixing 
angle parameters sin^ 612, sin^ 623 and sin^ ^13, the Dirac phase 6 and the rephasing 
invariant Jcp associated with 6, in the cases of the standard and inverse ordering 
of the charged lepton corrections. A summary and conclusions are presented in 
Section 5. Further details are reported in two appendices. In Appendix A we 
illustrate in detail the parametrisation we use for the standard ordering setup. 
Finally, in Appendix B we describe the statistical analysis used to obtain the 
numerical results. 

2 General Setup 

While neutrino masses and mixings may or may not look anarchical, the hierarchy 
of charged lepton masses suggests an ordered origin of lepton flavour. Given the 
wide spectrum of specific theoretical models, which essentially allows to account 
for any pattern of lepton masses and mixings, we would like to consider here the 
consequence for lepton mixing of simple, general assumptions on its origin. As we 
have indicated in the Introduction, we are interested in the possibility that the 613 
mixing angle originates because of the contribution of the charged lepton sector to 
lepton mixing [10, 11, 12, 13, 14, 15, 16]. The latter assumption needs a precise 
definition. In order to give it, let us recall that the PMNS mixing matrix is given 
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by 

U = UlU,, with ?7e, f/. defined by ; % 9 



V V 



where and rriy are respectively the charged lepton and neutrino Majorana mass 
matrices (in a basis assumed to be defined by the unknown physics accounting for 
their structure) and m^^^ and m^^^ are diagonal with positive eigenvalues. 

We will assume that the neutrino contribution Uy to the PMNS matrix U has 
[7^3 = 0, so that the PMNS angle ^13 vanishes in the limit in which the charged 
lepton contribution f/g can be neglected, f/g = 1- This is a prediction of a number 
of theoretical models. As a consequence, Vy can be parameterized as 

= ^uR23{e^23)Rl2{0u)^., (10) 

where Rij{0) is a rotation by an angle 6 in the ij block and "^y, $1, are diagonal 
matrices of phases. We will in particular consider specific values of d^^^ ^^d, in 
certain cases, of 62^, representing the predictions of well known models. 

The above assumption on the structure of Uy is not enough to draw conclusions 
on lepton mixing: any form of U can still be obtained by combining Uy with an 
appropriate charged lepton contribution Ue = U^W . However, the hierarchical 
structure of the charged lepton mass matrix allows to motivate a form of Ug similar 
to that of Uy, with U^^ = 0, so that we can write: ^ 

u, = ^,R2iie^23)Ri2ioi2)^e. (11) 

In fact, the diagonalisation of the charged lepton mass matrix gives rise to a 
value of that is small enough to be negligible for our purposes, unless the hier- 
archy of masses is a consequence of correlations among the entries of the charged 
lepton mass matrix or the value of the element {mE)3i, contrary to the common 
lore, happens to be sizable. In such a scheme, with no 13 rotation neither in the 
neutrino nor in the charged lepton sector, the PMNS angle ^13 is generated purely 
by the interplay of the 23 and 12 rotations in eqs. (10) and (11). 

While the assumption that Uf^ is small, leading to eq. (11), is well motivated, 
textures leading to a sizeable Ui^, are not excluded. In such cases, it is possible to 
obtain an "inverse ordering" of the R12 and -R23 rotations in Ue'. 

u, = ^,R^iiei2)R2im3)^e. (12) 

In the following, we will also consider such a possibility. 



^The use of the inverse in eqs. (11) and (12) is only a matter of convention. This choice 
allows us to lighten the notation in the subsequent expressions. 
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2.1 Standard Ordering 

Consider first the standard ordering in eq. (11). We can then combine U,^ and 
Ue in eqs. (10) and (11) to obtain the PMNS matrix. When doing that, the 
two 23 rotations, by the ^23 ^23 angles, can be combined into a single 23 
rotation by an angle 623- The latter angle is not necessarily simply given by the 
sum 6*23 = 6*23 + ^23 bccausc of the possible effect of the phases in \l/e (see 
further, eq. (68)). Nevertheless, the combination i?23(6'|3)^I^e^i/-R23(6'23) entering 
the PMNS matrix is surely a unitary matrix acting on the 23 block and, as such, 
it can be written as f2;^/223(^23)^e, where Q^^e are diagonal matrices of phases 
and 023 e [0,7r/2]. Moreover, we can write n,^i?23(^23)f^e = ^'u^R23{023)^'e, where 
$ = diag(l, e*"^, 1) and Q'^^ are diagonal matrices of phases that commute with 
the 12 transformations and either are unphysical or can be reabsorbed in other 
phases. The PMNS matrix can therefore be written as [15] 

U = PR,2{0l2)^R23{d23)Rl2{0l2)Q, (13) 

where the angle §23 can have any value, P is a diagonal matrix of unphysical phases, 
Q contains the two Majorana CPV phases, and $ = diag(l, e**^, 1) contains the 
only Dirac CPV phase. The explicit relation between the physical parameters 
6'23, and the original parameters of the model (^235 ^23' ^^"^ ^^'^ phases in 
\1/ = \l/*\l/!y) can be useful to connect our results to the predictions of specific 
theoretical models. We provide it in Appendix A. 

The observable angles in the standard PMNS parametrisation are given by 

sin 6*13 = \Ue3\ = sin 6112 sin ^23, 

• 2n \U^^3f . 2n COs"^ Of^ 

sm 6'23 = ; f-—2 = sm 623 7- 

1-\U,3\ 1 - sin^ 0^2 sin' ^23 (14) 

2 

I fj^ 1 2 sin e'(2 cos ei2 + e'^ cos O'^^ cos ^23 sin ^^2 
sin' 9i2 = = ^ . 

l-\Ue3\ 1 - sin' ^^2 sin' ^23 

The rephasing invariant related to the Dirac CPV phase, which determines the 
magnitude of CP violation effects in neutrino oscillations [21], has the following 
well known form in the standard parametrisation: 

Jcp = Im {U*^U*^Ue3U^i} = ^ sin 5 sin 26^13 sin 26^23 sin 26112 cos 6113 . (15) 

At the same time, in the parametrisation given in eq. (13), we get: 

Jcp = — sin sin 26'^2 sin 2^23 sin ^23 sin 2^5|'2 . (16) 
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sin ~023 
0.65 h 



0.60 : 
0.55 - 




Figure 1: The thick red line corresponds to the relation in eq.(17). The black 
and green lines show the la, 2a, 3a contours (solid, dashed and dotted lines, 
respectively) for sin ^13 and sin^ 623, as obtained in [8] (see Table 1). 



The relation between the phases and 6 present in the two parametrisations 
is obtained by equating eq. (15) and eq. (16) and taking also into account the 
corresponding formulae for the real part of U*^U*^Ue3U^i. To leading order in 
sin6'i3, one finds the approximate relation 5 ~ —(j) (see further eqs. (27), (28) and 
eqs. (32) and (33) for the exact relations). 

In this work we aim to go beyond the simplest cases considered already, e.g., in 
[15], where the charged lepton corrections to neutrino mixing are dominated only 
by the angle 6I2 and 623 is fixed at the maximal value 623 = 7r/4, and consider 
the case in which 623 is essentially free. A deviation of 623 from 7r/4 can occur in 
models in which ^23 = 7r/4 (BM, TBM) because of the charged lepton contribution 
to 623, or in models in which ^23 itself is not maximal (LC). This choice allows to 
account for a sizeable deviation of 623 from the value 7r/4, which appears to be 
suggested by the data [8]. If we keep the assumption ^^23 = 7r/4, the atmospheric 
mixing angle would be given by 

-j^ 2 gi]2^ 6 \ 1 

sin^ 6123 = -, — -0(1- sin^ ^13) , where sin 6113 = — = sin 6^12 • (17) 

2 1 - sm^ 6^13 2 ^J2 

This in turn would imply that the deviation from maximal atmospheric neutrino 
mixing corresponding to the observed value of ^13 is relatively small, as shown in 
Fig. 1. As for the neutrino angle 9^2-, we will consider two cases: 

TT 

• bimaximal mixing (BM): = ^ ^-Iso predicted by models with approx- 
imate conservation of L' = — — Lr ); 

• tri-hiraaxiraal mixing (TBM): = sin^^ —=. 



8 



Since in the approach we are following the four parameters of the PMNS matrix 
(the three measured angles 6*12, 6^23, 9i3 and the CPV Dirac phase S) will be ex- 
pressed in terms of only three parameters (the two angles 9I2, ^23 and the phase 
0), the values of 612, 623, O13 and 6 will be correlated. More specifically, 6 can be 
expressed as a function of the three angles, 6 = 5 (6^12, 6^23, 6*13), and its value will 
be determined by the values of the angles. As a consequence, the Jcp factor also 
will be a function of 612, 623 and ^13, which will allow us to obtain predictions for 
the magnitude of the CP violation effects in neutrino oscillations using the current 
data on sin^6'i2, sin^ 6*23 and sin^^ 



sin^^ 



13- 

We note first that using eq. (14) we can express sin^ ^23 in terms of sin^ 6'23 and 
13- 

.20 Sin^ ^23 - Sin^ ^13 2n COs2^23 ..gx 

sm ^23 = — , cos ^23 = • (18) 

1 — sm t7i3 1 — sm ^13 



It follows from these equations that ^23 differs little from 6'23 (it is somewhat larger). 
Further, using eqs. (14) and (18), we can express sin^ ^12 in terms of 9i2, O23, O13 
and 0: 

sin^ 6112 = (1 - cos^ 6^23 cos^ 6*13)"^ [sin^ e'^2 sin^ ^23 + cos^ 6'(2 cos^ 6^23 sin^ 6113 

+ 1 sin 20^2 sin 2^23 sin ^13 cos (p] . (19) 



As we have already indicated, we will use in the analysis which follows two specific 
values of e'(2 = 7r/4 (BM or LC); sin-i(l/y3) (TBM). Equation (19) will lead in 
each of the two cases to a new type of "sum rules", i.e., to a correlation between 
the value of 6*12 and the values of 623, O13 and (p. In the case of bimaximal and 
tri-bimaximal $'(2, the sum rules have the form: 

T3A/r -2/1 1 1 sin 2^23 sin ^13 cos 

BM : sm 612 = - + - ^ 20 

2 2 1 — COS^ 6'23 COS^ 6'i3 

= ^ + cot ^23 sin ^13 COS 0(1- cot^ ^23 sin^ + C»(cot^ ^23 sin^ ^13)) , 



(21) 



2 1 ( o sin 2^23 sin 6^13 cos — sin^ 6*23 



TBM: sin^^i2 = ;. 2 + ^^ ^ ^- ^ (22) 

3 \ 1 — COS^ 6^23 COS"' 0i3 ' 

= ]- 1 + 2V2 cot 6*23 sin 6*13 cos (l - cot^ 6*23 sin^ 6*13) 

+ cot^ ^23 sin^ + C(cot^ ^23 siu^ ^13)] . (23) 

The expressions for sin^ 612 in eqs. (20) and (22) are exact, while those given in eqs. 
(21) and (23) are obtained as expansions in the small parameter cot^ 6'23 sin^ 6*13. 
The latter satisfies cot^ 6'23 sin^ 6*13 < 0.063 if sin^ 6'23 and sin^ 613 are varied in the 
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3(7 intervals quoted in Table 1. To leading order in sin6'i3 the sum rule in eq. (21) 
was derived in [12]. 

We note next that since 6*12, 6^23 and 6*13 are known, eq. (19) allows us to express 
COS0 as a function of 612, O23 and 6is and to obtain the range of possible values of 
0. Indeed, it follows from eqs. (20) and (22) that 

RA/r A cos 2gi2(l-COS^ 023 cos^ ^13) 

BM : cos = — — , 24 

sm/t^23 smt^i3 

rr.^.. , (3sin2 0i2 -2) (1 -COS2 023COS2 013) + sin2 023 , x 

TBM : cos0=^ — -. (25) 

V2sin26'23 sin^^i3 

Taking for simplicity for the best fit values of the three angles in the PMNS matrix 
sin^6'i2 = 0.31, sin^ 6^23 = 0.39 and sin6'i3 = 0.16 (see Table 1), we get: 

COS0 ^ -0.99 (BM); cos ^ - 0.20 , (TBM). (26) 

Equating the imaginary and real parts of U*iU*^Ue3U^i in the standard para- 
metrisation and in the parametrisation under discussion one can obtain a relation 
between the CPV phases S and 0. We find for the BM case (^^^^2 ~ 7r/4): 

smd = — — , 27 

sin 2^12 ^ ^ 

cos 4> f 2 sin^ 6*23 

sin 2^ 

Since, as can be easily shown, 

/ cot^ 023 sin^ 6'i 3 cos^ \ ^ , , 

sin 2^12 = 1 - 4 . , , 29 

V (1 +COt2e23Sm^^i3)V 

we indeed have to leading order in sin6'i3, sin 5 = — sin0 and cos 5 = cos (p. 

The expressions for sin 5 and cos 5 in eqs. (27) and (28) are exact. It is not 
difficult to check that we have sm^ 6 + cos^ 5 = 1. Using the result for cos0, eq. 
(24), we can get expressions for sin 5 and cos 5 in terms of 6*12, 6*23 and 6*13. We 
give below the result for cos 6: 

cos5 = - / cot 26*12 tan 6*23 (l - cot^ 6*23 sin^ 6113) . (30) 
2 sin C7i3 

Numerically we find for sin^ 6*12 = 0.31, sin^ ^23 = 0.39 and sin 6*13 = 0.16: 

sin5 = ±0.170, cos5^ -0.985. (31) 



cos(5 = ^ — ^ ^ 1 . (2^ 

sin 2^12 Vsm^^23Cos2^i3 + sin2ei3 ' ^ 
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Therefore, we have 5 ~ vr. For fixed sin^ 6'i2 and sin 6*13, | cos5| increases with the 
increasing of sin^6'23. However, sin^ ^23 cannot increase arbitrarily since eq. (20) 
and the measured values of sin^ 9i2 and sin^ ^13 imply that the scheme with bi- 
maximal mixing under discussion can be self-consistent only for values of sin^ ^23, 
which do not exceed a certain maximal value. The latter is determined taking into 
account the uncertainties in the values of sin^ 9i2 and sin ^13 in Section 3, where 
we perform a statistical analysis using the data on sin^6'23, sin^6'i2, sin ^13 and 6 
as given in [8]. 

In a similar way we obtain for the TBM case {612 = sin^^(l/A/3)): 

. 2V2 sin0 

smS= ■ 00 ^ ^2 

3 sm26'i2 

2^2 , f 2sin2^23 

COSO = — - — — — COS 



3 sin 2^12 ^ " \^ sin2 ^^^2 q^^ _^ gjj^a 

1 sin 26^23 sin 6^13 

3 sin 2^12 sin^ ^23 cos2 ^13 + sin^ 6^^' ^ ' 

The results for sin 5 and cos 5 we have derived are again exact and, as can be 
shown, satisfy sin^ 5 + cos^ 5 = 1. Using the above expressions and the expression 
for sin^ 612 given in eq. (22) and neglecting the corrections due to sin 6*13, we obtain 
sin 5 ~ — sin0 and cos 5 ~ cos(^. With the help of eq. (25) we can express sin 5 
and cos 5 in terms of 612, O23 and 613. The result for cos 5 reads: 

^°^'- 3.^WM [l-.(3.n^...-2) (l-cot^fe.n^9.)] • (34) 
For the best fit values of sin^ = 0.31, sin^ 6*23 = 0.39 and sin 6*13 = 0.16, we find: 

sin5 = ±0.999, cos 5 ^ - 0.0490 . (35) 

Thus, in this case S ^ 7i/2 or 37r/2. For sin^ 6*23 = 0.50 and the same values of 
sin^^i2 and sin^ ^13 we get cos 5 = —0.096 and sin 5 = ±0.995. 

The fact that the value of the Dirac CPV phase S is determined (up to an 
ambiguity of the sign of sin 5) by the values of the three mixing angles 6*12, ^23 and 
6'i3 of the PMNS matrix, eqs. (30) and (34), are the most striking predictions of 
the scheme considered with standard ordering and bimaximal and tri-bimaximal 
mixing in the neutrino sector. For the best fit values of 612, O23 and 613 we get 
6 = TT and 5 = 7r/2 or 37r/2 in the cases of bimaximal and tri-bimaximal mixing, 
respectively. These results imply also that in the scheme with standard ordering 
under discussion, the Jcp factor which determines the magnitude of CP violation 
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in neutrino oscillations is also a function of the three angles 612, 623 and Oi^ of the 
PMNS matrix: 

JcP = Jcp{(^12, 6*235 6*13, S{6i2, O23, 9is)) = JcpiOu, O23, ^13) • (36) 

This allows to obtain predictions for the range of possible values of Jcp using 
the current data on sin^6'i2, sin^6'23 and sin ^13. We present these predictions in 
Section 3. The predictions we derive for 6 and Jcp will be tested in the experiments 
searching for CP violation in neutrino oscillations, which will provide information 
on the value of the Dirac phase S. 



2.2 Inverse Ordering 

As anticipated, we also study for completeness the case where the diagonalisation 
of the charged lepton mass matrix gives rise to the inverse ordering in eq. (12). 
The PMNS matrix, in this case, can be written as [11] 

U = i?23(^l3)^12(^12)^^23(^23)^12(^r2)0, (37) 

where unphysical phases have been eliminated, Q contains the two Majorana 
phases, and \E' = diag(l, e*'^, e*"^). Unlike in the case of standard ordering, it is 
not possible to combine the 23 rotation in the neutrino and charged lepton sec- 
tor and describe them with a single parameter, ^23- After fixing ^23 ^^id 612, we 
therefore have, in addition to the Majorana phases, four independent physical pa- 
rameters, two angles and two phases, one more with respect to the case of standard 
ordering. In particular, it is not possible anymore to write the mixing matrix in 
terms of one physical Dirac CPV phase only. Thus, in this case the four param- 
eters of the PMNS matrix (the three angles 612, O23 and 613 and the Dirac CPV 
phase 6) will be expressed in terms of the four parameters of the inverse ordering 
parametrisation of the PMNS matrix, eq. (37). We have for sin ^13, sin ^23 and 
sin 6 



12- 



sin 6^13 — SI2S23, 

• n V \V'23) ^23^^ '-12'-23 

•Smd23 = S23 — 



sin 6 



12 - 612- 



VI - Csush? ' (3^ 

e „u I 
12 '-23 



\pie _|_ piip(fu ^-lse ; 
\^12^ ^ \^12) ^12'^: 



1 (■^12-^23)^ 
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Given that the expressions for 6*23 and ^^13 do not depend on the value of 612, they 
will be the same for bimaximal and tri-bimaximal mixing (in both cases 6^^ = |): 

sin^,3 = ^, (39) 

2„ 1 1 + sin 26^20 v^cos"2^ COS — 2 sin^ 6*13 cos^ 6^20 . , 

sm 023 = ^ 40 

2 cos^ 6*13 

= ^ (^1 + sin 26^23 cos co' - cos 29^23 sin^ ^13 + C'(sin^ ^13) ) , (41) 

where the phase u' = ip — u. The expression (41) for sin^ 623 is approximate, the 
corrections being of the order of sin^ 613 or smaller. 

For each value of the phase ip, any value of ^13 and 623 in the experimentally 
allowed range at a given C.L., can be reproduced for an appropriate choice of u', 
^Ig. This is not always the case for the solar neutrino mixing angle ^12, as 
we will see in Section 4. Using eqs. (39), sin^ 6*12 can be expressed in terms of 613 
and if) as follows: 

vr 

• bimaximal mixing (BMio), = ~r- 



sin^ 9i2 = - — f 1 + 2 sin 613 v/cos2^ cos - sin^ ^13) (42) 

2 cos"^ 613 \ / 

~ - + sin 6^13 cos ip + C(sin^ 6^13) ; (43) 
2 

• tri-bimaximal mixing (TBMio), 9 12 = sin~^ 

sm^9i2 = - — fl + 2V2sin^i3v/cos20i3Cos^) (44) 

3 cos^ 0i3 V / 

~ ^ (1 + sin^ ^13) + ^ sin ^13 cos + C(sin^ ^13) . (45) 

The expressions for sin^ 9i2 in eqs. (42) and (44) are exact, while those given 
in (43) and (45) are obtained as expansions in sin^ ^13 in which the terms up to 
C(sin^6'i3) and C(sin^6'i3), respectively, were kept. Note that the corrections to 
the approximate expressions for sin^ 6*12 are negligibly small, being C(sin'^ 6*13). 
This together with eq. (43) and the 3a ranges of allowed values of sin^ 9i2 and 
sin6'i3 quoted in Table 1 suggests that the bimaximal mixing scheme considered 
by us can be compatible with the current (3cr) data on sin^ 9i2 and sin ^13 only for 
a very limited interval of negative values of cosV' close to (— 1). 
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It follows from eqs. (42) and (44) that the value of cosip is determined by the 
values of the PMNS angles 9i2 and ^13. At the same time, sin^ 6*23 depends on two 
parameters: u' and ^23- This implies that the values of u' and ^23 ^ire correlated, 
but cannot be fixed individually using the data on sin^ 6'23. 

It is not difficult to derive also the expressions for the Jcp factor in terms of 
the inverse ordering parameters in the two cases of values of 6'(2 of interest: 



BM : Jcp ~ ^ [sin cos 29^^ + sin u cos sin 2^^3 j + 0{sin^ ^13) , (46) 

TBM: Jcp ^ - ^(sin^ cos 26^23 + ^' cos ip sin 26^23 j + C'(sin^ 6^13) . (47) 

We have not discussed here the LC case (conservation of the lepton charge 
L' = Le — — Lr) as it involves five parameters (6'23, 6'f2, 6*235 ^.nd two CPV 
phases). At the same time, the "minimal" LC case with ^23 = is equivalent to 
the standard ordering case with BM mixing (i.e., with = '^/4) analised in detail 
in the previous subsection. 

As in the case of the standard ordering, to obtain the CPV phase 5 of the 
standard parametrisation of the PMNS matrix from the variables of these models, 
that is the function 5 = 6{ip, 6*23, 6*13), we equate the imaginary and real parts of 
^ei^^s^esUfj,! in the two parametrisations. 



3 Results with Standard Ordering 

In the numerical analysis presented here, we use the data on the neutrino mixing 
parameters obtained in the global fit of [8] to constrain the mixing parameters 
of the setup described in Section 2. Our goal is first of all to derive the allowed 
ranges for the Dirac phase 6, the Jcp factor and the atmospheric neutrino mixing 
angle parameter sin^ 6'23. We will also obtain the allowed values of sin^ 9i2 and 
sin^^i3. We start in this Section by considering the standard ordering setup, and 
in particular the two different choices for the angle 9'(2'- 9^2 = 7r/4 (BM and LC), 
9'(2 = sin-i(l/\/3) (TBM). 

We construct the likelihood function and the for both schemes of bimaximal 
and tri-bimaximal mixing as described in Appendix B, using as parameters for this 
model sin ^13, sin^ 6'23 and 6, and exploiting the constraints on sin^^i2, sin^6'23, 
sin^ ^13 and on 6 obtained in [8] . 

In Fig. 2 we show the contours of N„ = \f)^ in the (sin^6'23,5) plane, where 
the value of sin 6*13 has been marginalized. The blue dashed lines represent the 
contours of constant Jqp (in units of 10~^). In Figs. 3 and 4, starting from the 
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Standard Ordering - Normal Hierarchy 
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Figure 2: Contour plots for 




6 >6 



X in the standard ordering setup and 



normal hierarchy of neutrino masses. The value of the reactor angle ^13 has been 
marginalized. The solid, dashed and dotted thick lines represent respectively the 
1(7,2(7 and 3(7 contours. The dashed blue lines are contours of constant \ Jcp\ in 
units of 10"^. 



same likelihood function, we show the bounds on the neutrino mixing parameters 
and Jcp in each scheme, both for normal and inverted neutrino mass hierarchy. 
These bounds are obtained minimizing the in the parameter space of the model, 
keeping as a constraint the value of the corresponding parameter. To make a 
direct comparison of the bounds obtained in the scheme considered by us with the 
general bounds obtained in the global fit in [8], we show the results from [8] with 
thin dashed lines. Thus, the thin dashed lines in Fig. 4 are the bounds on Jcp 
obtained using directly the results of the global fit [8] and eq. (15), and represent 
the present status of our knowledge on this observable assuming the standard 3- 
neutrino mixing setup ^. The thick solid lines represent the results obtained in the 
scheme with standard ordering considered. The blue and red color lines correspond 
respectively to the cases of normal and inverted neutrino mass hierarchy; in the 

^More refined bounds on Jcp in the standard parametrisation of the PMNS matrix could be 
obtained by the authors of [8], using the full likelihood function. 
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Standard Ordering - Tribimaximal 




Standard Ordering - Bimaximal 




sin (?|3 sin"623 sin"6i2 



Figure 3: No- as a function of each mixing angle for the TBM and BM models 
in the standard ordering setup. The dashed lines represent the results of the 
global fit reported in [8] while the thick ones represent the results we obtain in 
our setup. Blue lines are for normal hierarchy while the red ones are for inverted 
hierarchy (we used purple when the two bounds are approximately identical). 
These bounds are obtained minimizing the value of N^- in the parameter space 
for fixed value of the showed mixing angle. 
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Standard Ordering - Tribimaximal 




6 -0.04 -0.02 0.00 0.02 0.04 



Standard Ordering - Bimaximal 




Figure 4: No- as a function of 5 and Jcp for the TBM and BM models in the 
standard ordering setup. The dashed Hues represent the results of the global fit 
reported in [8] while the thick ones represent the results we obtain in our setup. 
Blue lines are for normal hierarchy while the red ones are for inverted hierarchy. 
These bounds are obtained minimizing the value of in the parameter space 
for a fixed value of 6 (left plots) or Jcp (right plots). 
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case when the two bounds are essentially identical we used purple color lines. 

From Figs. 2 and 3 we see that both the tribimaximal and bimaximal cases are 
well compatible with data. The la difference between the minimum of N^r in the 
two cases is due to the fact that the bound on 6 obtained in [8] favours values of 
5 ~ vr (see Table 1), which is indeed the value needed in the bimaximal mixing 
(or LC) scheme to lower the value of 612 from ^5^2 = while the tri-bimaximal 
mixing scheme prefers | cos6\ <^ 1 (see Subsection 2.1). 

The results we obtain for sin^6'i2, sin^ 6'23 and sin^ 6'i3 (i.e., the best fit values 
and the 3a ranges) in the case of tri-bimaximal mixing are similar to those given 
in [8]. In contrast, our results for the Dirac phase 6 and, correspondingly, for the 
Jcp factor, are drastically different. For the best fit values and the 3a allowed 
ranges of 6 and Jcp we find (see also Table 2): 



NH: 5^4.64 = y, 1.38 <(5< 1.97, or (48) 

4.29<5<4.91, (49) 

371 

IB.: 6^ 4.64 = — , 1.39 <6< 2.17 , or (50) 

4.04 < (5 < 4.93 , (51) 

NH : Jcp = -0.034 , 0.028 < Jcp ^ 0.039 , or (52) 

-0.039 < Jcp ^ -0.028 , (53) 

IH : Jcp = -0.034 , 0.027 < Jcp ^ 0.039 , or (54) 

-0.039 < Jcp ^ -0.026 . (55) 



The 3a intervals of allowed values of S (Jcp) in eqs. (48) and (50) (eqs. (52) and 

(54) ) are associated with the local minimum at 5 = 7r/2 [Jcp — 0.034) in Fig. 4 
upper left (right) panel, while those given in eqs. (49) and (51) (eqs. (53) and 

(55) ) are related to the absolute minimum at 5 = 37r/2 (Jcp = —0.034). 

The results we have obtained, reported in Figs. 2 and 4, and in eqs. (48) - 
(55), are quasi-degenerate with respect to Jcp —Jcp, or 5 — )■ {2tt — 6). This 
stems from the fact that the phase enters into the expressions for the mixing 
angles only via its cosine, see eqs. (20) and (22). This symmetry is slightly broken 
only by the explicit bound on 6 given in [8], which is graphically represented in 
Fig. 4 by the asymmetry of the dashed lines showing that negative values of Jcp 
are slightly favored. 

As Figs. 2 and 4 show, in the case of tri-bimaximal mixing, the CP conserving 
values of 5 = 0; tt; 27r is excluded with respect to the best fit CP violating values 

■^These ranges are obtained imposing: ^/Ax^ — \/N^ — (Af™"i)2 ^ 3. 
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Best fit 


3(T range 




Jcp (NH) 


-0.034 


-0.039 -. 


- -0.028 © 0.028 ^ 0.039 




Jcp (IH) 


-0.034 


-0.039^ 


- -0.026 0.027^0.039 




5 (NH) 


4.64 


1.38 


-^ 1.97 4.29 ^4.91 


TBM 


5 (IH) 


4.64 


1.39 


-^ 2.17 ©4.04 ^4.93 




sin ^^13 


0.16 




0.13 0.18 




sin^ (NH) 


0.39 




0.33^0.64 




sin2 (IH) 


0.39 




0.34^0.66 




sin^ 9i2 


0.31 




0.25^0.36 




Jcp 


0.00 




-0.027^0.026 




5 (NH) 


3.20 




2.35^3.95 




5 (IH) 


3.27 




2.37^3.94 


BM 


sin 6*13 


0.16 




0.13^0.18 




sin2 ^23 (NH) 


0.38 




0.33^0.47 




sin2 ^23 (IH) 


0.39 




0.34^0.50 




sin^ 9i2 


0.31 




0.28^0.36 



Table 2: Best fit and So" ranges in the standard ordering setup. When not 
explicitly indicated otherwise, the result applies both for normal hierarchy and 
inverted hierarchy of neutrino masses. 



5 = 7r/2;37r/2 at more than Aa. Correspondingly, Jcp = is also excluded with 
respect to the best-fit values Jcp — (—0.034) and Jcp — 0.034 at more than ia. 
It follows from eqs. (48) - (55) (see also Table 2) that the 3a allowed ranges of 
values of both 6 and Jcp form rather narrow intervals. These are the most striking 
predictions of the scheme with standard ordering and tri-bimaximal mixing under 
investigation. 

We obtain different results assuming bimaximal mixing in the neutrino sector. 
Although in this case the best fit values of sin^ Ou, sin^ 6*23, sin^ 6*13 and 6 practically 
coincide with those found in [8] , the 3a allowed intervals of values of sin^ 612 and 
especially of sin^ ^23 and 6 differ significantly from those given in [8]. 

For the best fit values and the 3a intervals of sin^ 612 and sin^ ^23 we get (see 
also Table 2): 

sin^^is = 0.31, 0.28 < sin^ ^12 < 0.36; (56) 

NH : sin^ 623 = 0.38 , 0.33 < sin^ ^23 ^ 0.47 ; (57) 

IH : sin^ 623 = 0.39 , 0.34 < sin^ 623 ^ 0.50 . (58) 
As in [8], we find for the best fit value of 6 and Jcp'- S = n and Jcp = 0. However, 
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the 3(7 range of 6 and, correspondingly, of Jcp, we obtain differ from those derived 
in [8]: 



NH: 2.35 <5< 3.95; - 0.027 < J^p ^ 0.026 . (59) 
IH: 2.37<5<3.94; - 0.027 < J^p ^ 0.026 . (60) 

We see, in particular, that also in this case the Dirac CPV phase 6 is constrained 
to lie in a narrow interval around the value 5 ~ vr. This and the constraint 
sin^^23 ^ 1/2 are the most important predictions of the scheme with standard 
ordering and bimaximal neutrino mixing. 



4 Results with the Inverse Ordering 

The case of inverse ordering is qualitatively and quantitatively different from the 
case of standard ordering. For given values of 612, ^235 t^e number of parameters 
is the same as in the PMNS matrix. Still, not all values of U can be obtained, as 
we shall see. 

The constraints on the reactor and atmospheric neutrino mixing angles are the 
same for bimaximal and tri-bimaximal mixing and can be derived directly from eq. 
(39). For any given value of the phase ip, any values of 6*13 and ^23 in the ranges 



0<sin^i3< ^, 



1 



(61) 

COs4 



O<sin2 023<^^^=^l + 0(sm^^^i3 



can be obtained by an appropriate choice of u', and ^23- Clearly, the range of 
values allowed for ^13 and ^23 covers the full experimentally allowed range. The 
solar neutrino mixing angle can now be expressed in terms of 613 and ip as in 
eq. (39). Any value of 612 in the interval 

BMio : 

1 1 — 2 sin^^i3-\/cos"2^ — sin^ ^13 2n 11 + 2 sin 0i3a/cos 26'i3 — sin^ ^13 

2 - '''' - 2 ' 

(62) 

TBM/o : 

11 - 2^2 sin 013 Vcos 2^13 ^ . 2n / 11 + 2^2 sin 013 Vcos 2^13 

Q ^ <sm0i2<- ^ , 63 

3 COS"' 013 3 COS^ 013 
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can then be obtained for an appropriate choice of ip. At leading order in sin^is 
these bounds become 



BMjo : ^ - sin ^13 < sin^ ^12 < ^ + sin ^13, 



(64) 

„ , , „ „ , , , - sin Bti. 

3 3 ~ ~ 3 3 



1 2 /2 1 2 /2 

TBM/o : - - ^ sin ^13 < sin^ < + ^ sin ^13. 



Given the experimental bounds on the PMNS angles found in the global fit [8], see 
Table 1, one can immediately notice that while the tri-bimaximal case is perfectly 
compatible with the data, the bimaximal case has a ~ 2cr tension in the prediction 
of the solar neutrino mixing angle parameter sin^ Ou- 

As was done for the standard ordering case, we construct the likelihood function 
and the for both models as described in Appendix B, exploiting the constraints 
on sin^^i2, sin^^23, sin^ ^13 and on 6 obtained in [8], and using in this case as 
parameters sin ^13, sin ^23 the phases ip and u. We show in Figs. 5 and 6 the 
bounds on the neutrino mixing angles and the Jcp factor both in the cases of 
bimaximal and tri-bimaximal mixing in the neutrino sector, and for normal and 
inverted neutrino mass hierarchy. 

From Fig. 5, we see that in the case of tribimaximal mixing (upper row), the 
intervals of allowed values of the PMNS mixing angles obtained in the model under 
discussion and in the global fit performed in [8] coincide. This is a consequence of 
the fact that the 4D parameter space of the model considered completely overlaps 
with the experimentally allowed parameter space in the PMNS parametrisation 
and therefore it does not give any additional constraint. It is consistent with the 
analytic bounds reported above as well. 

In the case of bimaximal mixing instead (Fig. 5 lower row), only a portion of 
the relevant PMNS parameter space is reachable, a fact that is reflected in the 
bounds on sin^6'i2 given in eq. (64). Values of 6*12 in the upper part of its present 
experimental range are favoured in this case. 

In both cases of tri-bimaximal and bimaximal mixing from the neutrino sec- 
tor, the bounds on sin^ 613 and sin^ 612 corresponding to the normal and inverted 
neutrino mass hierarchy are approximately identical, while they differ for the at- 
mospheric neutrino mixing angle and for the Jcp factor. 

Considering the expressions for Jcp in eqs. (46) and (47) and Fig. 6, we see 
that within ~ la from the best-fit point, every value in the ranges 

\JS^\<"^-0M, |J^n<^^- 0.038, (65) 

is allowed, where we have used the la upper bound on sin 6*13 from Table 1. As a 
consequence, we cannot make more specific predictions about the CP violation due 
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Inverse Ordering - Tribimaximal 




sin 613 sin-023 sin 612 



Inverse Ordering - Biniaximal 




sin 613 sin'^J23 sin-612 



Figure 5: A'^,- as a function of each mixing angle for the TBM and BM models 
with the inverse ordering setup. The dashed lines represent the results of the 
global fit reported in [8] while the thick ones represent the results we obtain in 
our setup. Blue lines are for normal hierarchy while the red ones are for inverted 
hierarchy (we use purple when the two bounds are approximately identical). 
These bounds are obtained minimizing the value of N^- in the parameter space 
for fixed value of the showed mixing angle. 
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Inverse Ordering 




Bimaximal 



-0.04 - 0.02 0.00 0.02 0.04 
(a) 



Jcp 




-0.04-0.02 0.00 0.02 0.04 



Jcp 



(b) 



Figure 6: N„ as a function of Jcp for tlie TBM and BM models in the inverse 
ordering setup. The dashed lines represent the results of the global fit reported in 
[8] while the thick ones represent the results we obtain in our setup. Blue lines are 
for normal neutrino mass hierarchy while the red ones are for inverted hierarchy. 
These bounds are obtained minimizing the value of in the parameter space 
for a fixed value of Jcp- 



to the Dirac phases 6 in this case. This is an important difference with respect to 
the standard ordering scheme where, in the tri-bimaximal mixing case, relatively 
large values of the \Jc'p\ factor lying in a narrow interval are predicted at 3a and, 
in the bimaximal mixing case, 6 is predicted to lie at 3a in a narrow interval around 
the value of 5 ~ tt. 
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5 Summary and Conclusions 



In this paper we considered the possibihty that the neutrino mixing angle 613 
arises from the interplay of 12 and 23 rotations in the neutrino {Uu) and charged 
lepton (f/g) contributions to the PMNS neutrino mixing matrix {U = UlU^)- We 
generalized previous work [15] in two directions. First, we considered two possible 
orderings of 12 and 23 rotations in [/g, the "standard", Ue ~ R23R12, and the 
"inverse" , Ue ~ RI2R23, while keeping the standard ordering in the neutrino sector, 
Uu ~ -^23-^12- Second, in order to be able to accommodate a possible deviation 
of the atmospheric neutrino mixing angle 623 from 7r/4, we allowed the charged 
lepton 23 rotation angle (and possibly the neutrino one, in the standard case) to 
assume arbitrary values. We considered the cases in which Uiy is in the bimaximal 
or tri-bimaximal form, or in the form resulting from the conservation of the lepton 
charge — — Lr (LC). We took, of course, all relevant physical CP violation 
(CPV) phases into account. 

The case of normal ordering turns out to be particularly interesting. The 
PMNS matrix can be parameterized in terms of the charged lepton and neutrino 
12 rotation angles, 9\2 and 9^2-, an effective 23 rotation angle, 623 ~ ^23, and a 
CPV phase 0. Once 6^2 is fixed to the bimaximal (LC) or tri-bimaximal value, the 
number of parameters reduces to three, and the Dirac phase 5 in the PMNS matrix 
can be predicted in terms of the PMNS solar, atmospheric and reactor neutrino 
mixing angles ^12, ^23 and 613. Moreover, the range of possible values of the PMNS 
angles turns out to be constrained. 

In the tri-bimaximal case, the Dirac CPV phase 5 is predicted to have a value 
5 ~ 7r/2 or 5 ~ 37r/2, implying nearly maximal CP violation in neutrino oscilla- 
tions, while in the bimaximal (and LC) case we find 5 ^ n and, consequently, the 
CP violation effects in neutrino oscillations are predicted to be small. The present 
data have a mild preference for the latter option (see Table 1 and, e.g.. Fig. 4). 
Moreover, 6'23 is predicted to be below 7r/4 in the bimaximal case, which is also in 
agreement with the indications from the current global neutrino oscillations data. 
In the set-up considered by us, the 6^23 > 7r/4 solution of the global fit analysis 
in [8] is disfavored. 

The case of inverse ordering is qualitatively and quantitatively very different. 
Fixing Uy to the bimaximal or tri-bimaximal form is not sufficient to obtain a 
prediction: the number of free physical parameters in this case is four - two angles 
and two CPV phases. Still, not all values of the four physical parameters in the 
PMNS matrix, 6'i2, 923-, 9i3 and 5, can be reached in this parameterization. In the 
tri-bimaximal case, the ranges of parameters that can be reached overlaps with 
the experimental ranges, so that no predictions can be made. In the bimaximal 
case, however, this is not the case. One obtains, in fact, the approximate relation 
sin^ 612 ^ 1/2— sin 613, which is barely compatible with the data. As a consequence. 



24 



i) there is a tension in the above relation that worsen the quahty of the fit, and 

ii) values of 612 in the upper part of its present experimental range are preferred. 
In both cases, no predictions for the Dirac CPV phase 6 can be made. We did 
not consider here the LC it involves, in general, five parameters, while its 
"minimal" version, corresponding to setting 623 = 0, is equivalent to the standard 
ordering case with BM mixing (i.e., with ^J'g = vr/4). 

The fact that the value of the Dirac CPV phase 6 is determined (up to an 
ambiguity of the sign of sin 5) by the values of the three PMNS mixing angles, 612, 
023 and 9i3, eqs. (30) and (34), are the most striking predictions of the scheme 
considered with standard ordering and bimaximal (LC) and tri-bimaximal mixing 
in the neutrino sector. As we have already indicated, for the best fit values of 
6'i2 ; ^23 and 6'i3 we get 5 = vr and 5 = 7r/2 or 37r/2 in the cases of bimaximal and 
tri-bimaximal mixing, respectively. These results imply also that in the scheme 
with standard ordering we have discussed, the Jcp factor which determines the 
magnitude of CP violation in neutrino oscillations, is also a function of the three 
mixing angles: Jcp = Jcp(6'i2, 6*23, 6*13, (^(^12, 6*23, ^'13)) = Jcp(6'i2, 6*23, ^'13). This 
allowed us to obtain predictions for the range of possible values of Jcp using the 
current data on sin^ 9i2, sin^ 6'23 and sin6'i3, which are given in eqs. (48) - (54) and 
eqs. (59) - (60). 

The predictions for sin^ 6'23, and for S and Jcp we have obtained in the scheme 
with standard ordering and bimaximal (or LC) or tri-bimaximal form of U^, will be 
tested by the neutrino oscillation experiments able to determine whether sin^ 623 ^ 
0.5 or sin^ 623 > 0.5, and in the experiments searching for CP violation in neutrino 
oscillations. 
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Appendix 



A Parametrisation of the PMNS matrix 

In the present Appendix we show how the parametrisation of eq. (13) follows from 
the ones in eqs. (10) and (11). We start by writing explicitly the PMNS matrix as 

u = ^iRi2{ei^)R2Mz)'^R2Mz)Ri2{e\^)^,, (66) 

where ^ = diag(l, e*'^, e*"^), without loss of generality. Any 2x2 unitary matrix 
V can be recast in the form V = PR{6)Q, where P = diag(e*'^i, e*'^^), Q = 
diag(l, e*'^^) and R{9) is a 2 x 2 rotation. We use this to write 

i?23(^l3)*^23(^23) = $'^^23(^23)^^, (67) 

where -R23(^23) is an orthogonal rotation in the 23 block with 

sin ^23 = I cos sin + e'^"'"^) sin ^^3 cos ^^3 1 , (68) 

$' = diag(e*'^i,e*'^2^e*'^3), and Q = diag(l, e*'^^^ e*'^^). An explicit solution for the 
angles in terms of the original parameters is 

01 = 0, (f)2 = 5c + 5s + ^-UJ, (j)3 = 0, 

1^2 = -(^s + 1^, 1^3 = -^c + 1^, 

where 

Ss = Arg (cos ^2^3 sin 9^,, + e'^^-^^ sin 9^,, cos 9^,,) , 
4 = Arg (cos ^^3 cos 9^3 - e^^-^-^) sin ^^3 sin 9^^) . 
Considering now also the Ri2{9i2) rotation, we obtain 

i?23(^23)^^i?12(^^r2) = <f"i?23(^23)i?12(^r2) Q\ (71) 

with $" = diag(l,e*'^^e^'^2) ^nd Q" = diag(l, 1, e^('^3-'^2))_ phases in Q" add 
to the ones in Q' and are Majorana phases. The ones in instead, add to the 
ones in 

= e''^' diag(l, e'^'''^-'^'^'^^\ ei(</'3-<^i+^2))_ ^^2) 

The phase in the 33 position commutes with Ri2{9\2)- Together with the overall 
phase 01, it will describe the unphysical phase matrix P in eq. (13): 

P = e'^^ diag(l, 1, e'^'^'^-^'^'^^^) . (73) 

We see that the only physical Dirac CP violating phase in this parametrisation is 
contained in the matrix $ = diag(l, e*"^, 1), with 

= 02 - 01 + W2 = + 4- (74) 
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Figure 7: Icr, 2(t, 3(t contours (respectively solid, dashed and dotted lines) of our 
global likelihood function in the (sin^ ^23; sin^ ^13) plane (left) and (sin^0i3,(^) 
plane (right), using the data for NH. These plots can be compared with Fig.l 
and Fig. 2 of [8] for NH. Undisplayed variables have been marginalized. 



B Statistical analysis 

In this appendix we describe the simplified statistical analysis performed to obtain 
the results. Our aim is to use the results of the global fit performed in [8] to 
assess how well each of the models introduced in the previous section can fit the 
data. In particular, we use the constraints on the PMNS angles ^13, ^12, ^23 and on 
the phase 8 for the normal hierarchy (NH) and inverted hierarchy (IH) cases, as 
derived in [8] . There, the results are reported by plotting the value of iV^ = \/ Ax^ 
(with Ax^ = id — Xmin) as a function of each observable, with the remaining ones 
marginalized away. We construct an approximate global likelihood from these 
functions as 

L,(a,) = exp (-^^1^) , m = nL,(a,), (75) 

where a = {sin^ Oi^, sin^ 6'23, sin^ 612, 6} are the observables relevant for our analy- 
sis, and we define 

x'{d) = -2logLid) (76) 

and Nfj[d) = a/ X^i^)- In using this procedure we loose any information about 
possible correlations between different observables. The effect of this loss of in- 
formation is however negligible, as one can check comparing our la, 2a and 3a 
contours in the (sin^ ^23, sin^ ^13) and (sin^ 6*13, 5) planes shown in Figure 7 with 
the ones in Fig. 1 and Fig. 2 of [8]. 

Each model introduced in the previous section (which we dub with an index m) 
depends on a set of parameters x™ = {xf^}, which are related to the observables 
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via expressions aj = a^{x."^), obtained from eqs. (14), (38). We then construct 
the hkehhood function in the space of the parameters x™ as 

L'™(x'") = L(a'"(x™)). (77) 

We define x^(x'") = -21ogL"(x") and A/'<,(x™) = ^/^(x™). The last one is the 
function we use to produce the plots shown in Figures 2-6. Finally, to obtain the 
best-fit point we use the maximum likelihood method. 
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